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HIGHLIGHTS 


•  A  methodology  to  determine  the  modulus  of  functional  layers  in  SOFCs  is  reported. 

•  In-situ  Young's  modulus  for  various  functional  layers  is  reported. 

•  The  proposed  testing  methodology  could  be  applied  for  other  multilayer  systems. 
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In-situ  Young's  moduli  of  thin  constituent  layers  of  a  solid  oxide  fuel  cell  (SOFC)  are  needed  to  estimate 
the  mechanical  reliability  of  the  fuel  cell  unit.  Because  a  robust  technique  to  measure  the  same  is  not 
available,  an  improved  methodology  is  proposed  to  determine  the  in-situ  Young's  moduli  of  thin  ceramic 
layers  of  a  substrate-supported  SOFC.  The  measured  Young's  moduli  of  the  constituent  layers  were  found 
to  be  close  to  those  of  corresponding  bulk  materials  using  the  resonant  ultrasound  spectroscopy  (RUS) 
technique  but  were  different  from  the  values  obtained  using  nanoindentation. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  are  environmentally  clean,  high  ef¬ 
ficiency  functional  devices  that  transform  chemical  energy  stored  in 
fuel  into  electrical  energy  1-4  .An  SOFC  can  be  structurally  viewed 
as  a  multilayer  ceramic  composite  system  that  operates  at  elevated 
temperature.  The  individual  ceramic  layers  in  the  SOFC  contain  a 
multitude  of  microstructural  features  (e.g.  pores,  microcracks,  phase 
boundaries  etc.),  which  play  an  important  role  in  the  functional, 
transport  5,6]  and  mechanical  response  [7-11  of  the  cell.  For 
example,  an  increase  in  porosity  in  an  electrochemically  active  layer 
could  lead  to  a  voltage  drop  [  12  ,  decrease  in  electrolyte  temperature 
[5]  and  activation  overpotential  13]  and  may  lead  to  structural 
failure  as  a  result  of  the  degradation  of  Young's  moduli  of  thin 
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constituent  layers  5,14].  The  Young's  modulus  is  used  to  calculate  the 
hotspot  failure  stresses  that  are  generated  as  a  result  of  thermal 
gradients  across  the  ceramic  layer  [15-18  .  The  Young's  modulus  is 
also  required  to  determine  the  residual  stresses  that  develop  in  the 
individual  layers  as  the  SOFC  is  cooled  from  the  processing  temper¬ 
ature  and  subsequently  heated  back  to  the  operation  temperature 
[15,20-23  ].  A  robust  technique  to  measure  the  Young's  moduli  of  the 
constituent  layers  will,  therefore,  be  beneficial  for  the  mechanical 
design,  reliability  and  safe  operation  of  SOFCs. 

The  nanoindentation  technique  is  most  frequently  used  to 
obtain  the  Young's  modulus  of  thin  ceramic  materials  19,24,25]. 
For  specimens  with  high  defect  densities,  however,  the  material 
properties  obtained  from  the  nanoindentation  technique  are  often 
unreliable  because  they  constitute  the  local  instead  of  the  bulk 
response  of  the  material  [26  .  At  the  same  time,  in-situ  elastic 
properties  of  multilayer  systems  such  as  SOFCs  and  thermal  barrier 
coatings  (TBCs)  [27  have  proven  to  be  difficult  to  measure  exper¬ 
imentally,  owing  to  the  reduced  dimensionality  of  the  individual 
layers.  In  the  present  investigation,  an  improved  experimental 
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Fig.  1.  Schematic  of  substrate  supported  solid  oxide  fuel  cell  (SOFC).  The  substrate  is  a  rectangular  thick  wall  (<1  mm)  MMA  (magnesia-magnesium  aluminate  spinel)  material  that 
is  extruded  with  channels  for  fuel  flow.  Starting  from  the  porous  anode  barrier  (PAB),  subsequent  functional  layers  were  screen  printed  sequentially  onto  the  MMA  substrate. 


methodology  for  determination  of  the  Young's  modulus  of  multi¬ 
layer  materials  is  proposed.  The  current  methodology  consists  of. 

1.  Fabrication  of  a  series  of  specimens  with  increasing  number  of 
constituent  layers  that  are  screen-printed  on  a  porous  ceramic 
substrate. 

2.  Microtensile  testing  of  individual  bi-  and  multilayer  specimens 
under  ambient  conditions. 

3.  Measurement  of  in-plane  strain  on  both  sides  (substrate  and 
outermost  functional  layer  side)  of  the  specimen  during  tensile 
deformation  by  application  of  a  non-contact  digital  image  cor¬ 
relation  (DIC)  technique. 

4.  Calculation  of  the  Young's  modulus  of  each  thin  ceramic  con¬ 
stituent  layer  in  a  substrate  supported  SOFC  design  by  applying 
the  equistrain  criterion. 

The  above  approach  also  allowed  identification  of  the  validity 
limit  of  the  equistrain  criterion.  The  Young's  modulus  values  of  in¬ 
dividual  thin  ceramic  layers  determined  with  the  proposed  meth¬ 
odology  are  compared  with  measurements  performed  on  similar 
materials  with  nanoindentation  and  with  bulk  specimens  using  the 
resonant  ultrasound  spectroscopy  (RUS)  technique. 

2.  Experimental  methodology 

2  A.  Specimen  preparation:  multilayer  ceramic  specimens 

SOFC  samples  were  manufactured  by  LG  Fuel  Cell  Systems  Inc. 
(LGFCS).  Unlike  an  anode  supported  SOFC  [28  ,  the  structural  support, 
in  this  case,  is  provided  by  a  rectangular  thick  wall  (<1  mm)  MMA 
(magnesia-magnesium  aluminate  spinel)  extruded  substrate  with 


longitudinal  channels  for  fuel  flow  (Fig.  1).  Functional  layers  were 
screen-printed  in  a  sequential  manner  onto  the  MMA  substrate, 
starting  from  the  porous  anode  barrier  (PAB).  Additional  functional 
layers  that  were  screen  printed  on  the  substrate  in  sequence  included 
the  anode  current  collector  (ACC),  anode  (A),  electrolyte  (E),  cathode 
(C),  and  cathode  current  collector  (CCC).  This  sequential  screen 
printing  led  to  six  different  types  of  multilayer  samples;  each  with 
different  number  of  functional  layers  ( Table  1 ).  Microtensile  speci¬ 
mens  were  fabricated  from  the  substrate  and  multilayer  samples 
using  a  series  of  dry  cutting  and  grinding  procedures  described 
elsewhere  [29,30  .  Some  multilayer  specimens  had  a  curvature  prior 
to  cutting  and  the  grinding  process  reduced  the  specimen  curvature, 
presumably  due  to  residual  stress  relief  as  a  result  of  the  grinding 
operation.  In  these  multilayer  samples,  grinding  was  always  per¬ 
formed  on  the  substrate  side,  thereby  preserving  the  thickness  and 
integrity  of  the  functional  layers.  Further,  in  the  multilayer  samples, 
the  microtensile  specimens  were  fabricated  from  the  active  area 
of  the  fuel  cell  such  that  the  layers  were  all  continuous  in  the  plane  of 
the  specimens.  Microtensile  beam  specimens  were  prepared  to  a 
rectangular  test  geometry  with  a  total  average  thickness  of 
0.190  ±  0.002  mm,  gage  length  of  6  ±  2  mm  and  gage  width  of 
2  ±  0.01  mm  (  ^able  1 ).  In  addition,  MMA  substrate  specimens  with 
various  thickness  values  (190, 150,  and  100  pm)  were  prepared  for 
microtensile  testing  and  dimensions  of  the  specimens  were  measured 
optically  prior  to  the  testing  using  a  Nikon  profile  projector  (V-12). 

2.2.  Microstructure  and  testing  procedure 

The  microstructure  of  these  materials  was  characterized  using  a 
Hitachi  S-4800  Field  Emission  scanning  electron  microscope  (SEM). 
Uniaxial  tensile  experiments  were  conducted  on  the  MMA 


Table  1 

List  of  screen-printed  multilayer  specimen  types  and  their  composite  moduli. 


Multilayer  ceramic  specimen 

No  of  layers 

(excluding 

substrate) 

In-plane 

layers 

Net  thickness  of  substrate  + 
average  thickness  of  active  layers 
(total  thickness  of  ceramic 
composite  layers  -190  ±  5  pm) 

Figure  no# 

Composite  modulus 
(GPa)  using  20  pe 
deviation  from  linearity 

Substrate  only  (S) 

0 

S 

190,  170,  150,  130,  100,  50 

4 

70  ±  1 

Substrate  (S)  +  porous  anode  barrier  (PAB) 

1 

S  &  PAB 

165  +  25 

4a 

69  ±  2 

Substrate  (S)  +  porous  anode  barrier 
(PAB)  +  anode  current  collector  (ACC) 

2 

S  &  ACC 

150  +  25+15 

4b 

71  ±4 

Substrate  (S)  +  porous  anode  barrier 
(PAB)  +  anode  current  collector 
(ACC)  +  anode  (A) 

3 

S  &  A 

143  +  25+15  +  7 

4c 

71  ±  3 

Substrate  (S)  +  porous  anode  barrier 
(PAB)  +  anode  current  collector 
(ACC)  +  anode  (A)  +  electrolyte  (E) 

4 

S  &  E 

133  +  25+15  +  7+10 

4d 

79  ±  2 

Substrate  (S)  +  porous  anode  barrier 
(PAB)  +  anode  current  collector 
(ACC)  +  anode  (A)  +  electrolyte  (E)  +  cathode 
(C) 

5 

S  &  C 

113  +  25+15  +  7+10  +  20 

4e 

84  ±3 

Substrate  (S)  +  porous  anode  barrier 
(PAB)  +  anode  current  collector 
(ACC)  +  anode  (A)  +  electrolyte  (E)  +  cathode 
(C)  +  cathode  current  collector  (CCC) 

6 

S  &  CCC 

43  +  25+15  +  7+10  +  20+  70 

4f 

105  ±  2 

524 
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Fig.  2.  Cross-sectional  SEM  micrograph  of  the  substrate  supported  solid  oxide  fuel  cell 
(SOFC). 

substrate  and  on  all  multilayer  ceramic  specimens  using  a  micro¬ 
testing  system  described  elsewhere  [7].  Specimens  were  glued  at 
both  ends  of  the  tensile  fixture  and  were  deformed  at  a  displace¬ 
ment  rate  of  10-3  mm  s_1.  Images  of  the  gage  section  were  captured 
at  a  rate  of  s-1  using  a  high-resolution  CMOS  camera  and  strains 
were  calculated  with  the  digital  image  correlation  (DIC)  technique 
using  a  Matlab®-based  code  described  elsewhere  7,31,32  .  A  grid  of 
approximately  25,000  virtual  markers  was  defined  and  each 
marker  subset  was  20  x  20  pixels2.  Due  to  the  rectangular  geom¬ 
etry  associated  with  the  gage  section  of  the  microtensile  specimens, 
some  specimens  failed  at  the  glued  region-gage  section  interface. 

Resonant  ultrasound  spectroscopy  (RUS)  was  applied  to  mea¬ 
sure  the  elastic  modulus  of  disk  shaped  MMA  specimens 
(thickness  ~  0.70  mm  and  diameter  ~  31.8  mm).  RUS  was  also 
performed  on  anode,  electrolyte,  cathode  and  cathode  current 
collector  with  specimens  that  were  powder  pressed,  sintered  for 
1  h  at  1380  °C  and  then  machined  into  disks  (thickness  ~  1.5  mm 
and  diameter  ~  25  mm).  Nanoindentation  experiments  were  per¬ 
formed  to  an  indentation  depth  of  250  nm  using  a  Hysitron  Tri- 
boindentor™  with  a  Berkovich  tip  on  the  ceramographically 
prepared  cross-section  of  the  SOFC  sample.  A  total  of  15  in¬ 
dentations  were  made  in  each  SOFC  layer,  and  average  values  of 
reduced  modulus  [33]  are  reported. 

3.  Results  and  discussion 

3.1.  Microstructure 

Fig.  2  shows  an  SEM  micrograph  of  the  cross-section  of  a 
substrate  supported  SOFC  design  that  was  imaged  using 


3  Strain  (%) 


backscattered  electrons.  The  functional  layers  along  with  their 
individual  average  thickness  values  (in  the  range  7-70  pm)  are 
identified  in  Fig.  2.  Each  layer  in  Fig.  2  exhibited  varying  levels  of 
porosity,  with  the  electrolyte  exhibiting  the  lowest  level  of 
porosity.  The  microstructure  of  the  substrate  consists  of  a  porous 
composite  of  cubic  magnesium  oxide  (MgO)  and  spinel 
(MgAl204)  phases.  In  addition,  the  substrate  also  contained  an 
appreciable  density  of  engineered  defects  (microcracks  and 
porosity)  that  facilitate  the  functional  requirement  of  fuel  flow 
across  the  thickness  of  the  substrate.  Microcracking  in  the  MMA 
substrate  is  a  result  of  thermal  strain  mismatch  between  the  two 
cubic  phases  during  cooling  from  the  processing  temperature  and 
the  microcracks  are,  therefore,  referred  to  as  thermally  induced 
microcracks  [7  .  The  cell  structure  consists  of  multiple  layers:  the 
YSZ  (yttria-stabilized  zirconia),  porous  anode  barrier  (PAB),  a 
cermet  anode  current  collector,  the  active  anode,  a  ScSZ  elec¬ 
trolyte,  LSM  (lanthanum  strontium  manganate)-YSZ  based  cath¬ 
ode  and  an  LSM-based  cathode  current  collector.  It  is  noted  that 
the  unit  cell  structure  of  the  SOFC  has  a  functionally  graded 
structure  with  varying  thickness,  microstructure  and  density. 
This  design  is  engineered  to  compensate  the  thermal  and  me¬ 
chanical  mismatches  arising  from  the  differences  in  thermal 
expansion  and  elastic  modulus  between  the  adjacent  layers  while 
achieving  the  required  electrochemical  performance. 

3.2.  Uniaxial  tensile  response:  MMA  substrate  specimens 

Two  sets  of  uniaxial  tensile  experiments  were  conducted, 
which  included  monotonic  loading  (190  pm  thick,  Fig.  3a)  and 
loading-unloading-reloading  (for  190  and  100  pm  thick  speci¬ 
mens,  Fig.  3b).  These  experiments  were  conducted  to  study  the 
effect  of  specimen  thickness  on  the  stress-strain  response  and  to 
decouple  anelastic  (recoverable)  and  inelastic  (irrecoverable) 
strains  at  various  stress  levels.  The  stress-strain  response  of  the 
specimen  can  be  a  function  of  the  specimen  thickness  due  to  the 
fact  that  the  employed  grinding  procedure  affects  the  response 
of  thinner  specimens  more  than  the  thicker  specimens.  Obser¬ 
vations  from  such  load-unload-reload  experiments  are  used  to 
determine  the  stress  value  beyond  which  pre-existing  micro¬ 
cracks  start  to  propagate  [7,10,11  .  Fig.  3a  shows  the  uniaxial 
tensile  response  of  MMA  substrate  (190  pm  thick)  from  three 
separate  tests  and  the  average  failure  stress  and  strain  were 
9.5  ±  0.2  MPa  and  0.032  ±  0.005%,  respectively.  Unlike  con¬ 
ventional  structural  ceramics,  tensile  response  of  the  MMA 
substrate  is  non-linear,  which  implies  that  the  apparent  modulus 
becomes  a  function  of  stress  or  strain  [7  .  This  non-linear  tensile 
response  of  porous  microcracked  ceramic  materials  is  dictated 
by  the  pre-existing  microcracks  and  their  evolution  during 
loading  [7]. 


10 


Fig.  3.  a:  Monotonic  uniaxial  tensile  stress-strain  response  of  MMA  substrate  specimen,  b:  Loading-unloading  and  reloading  response  of  MMA  substrate  specimens  of  different 
thickness  under  uniaxial  tensile  loading. 
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Existing  experimental  techniques  used  to  study  the  elastic  properties  of  multilayer  structures. 


525 


Layers/materials 

Techniques 

Reference# 

Thermal  barrier  coating  (bilayer  beam, 
top  coat/bond  coat) 

Micro  flexure  testing  and  composite 
beam  theory 

Eberl  et  al.,  2011  [27] 

Bi  and  tri  layer  ceramic  structures 

Flexural  testing  and  composite  beam 
theory 

Borba  et  al.,  2011  [34] 

Thin  film  and  coatings 

Uniaxial  Tension 

Pandey  et  al.,  2013  [32],  Deng  et  al.,  [36],  Hollman  et  al.,  1997  36] 

YSZ/NiCoCrAlY  composite  coatings 

Four-point  bending,  composite  beam 
theory,  rule  of  mixture 

Khor  et  al.,  2001  [38] 

Multilayer  thin  films 

Electrostatic  pull-in  test  structure  and 
analytical  modeling 

Nie  et  al.,  2005  [39] 

Cr/CrN  multilayer  coatings 

Microindentation 

Kot  et  al.,  2008  [40] 

AlN/TiN  thin  films 

Nanoindentation 

Rawdanowicz  et  al.,  1999  [41] 

Thin  flexible  multi-layered  material 

Compression  column  method 

Ohtsuki,  2011  [42] 

Thin  film  and  coatings 

Diffraction  and  curvature  method 

Chang  et  al.,  2009  [43],  Chen  et  al.,  2004  44],  Matejicek  and  Sampath, 
2003  [45],  Tran  et  al.,  1997  [46] 

Bi  and  tri  layer  ceramic  structures 

Combination  of  FE  and  RUS  technique 

Araki  et  al.,  2013  [47] 

The  loading-unloading-reloading  response  for  190  and  100  pm 
thick  specimens  of  MMA  substrate  are  shown  in  Fig.  3b.  Unloading 
was  conducted  from  five  different  stress  levels  (2.5,  5,  7,  8  and 
9  MPa)  and  the  corresponding  residual  strain  was  obtained  at  zero 
load.  It  was  found  that  a)  the  unloading  tensile  response  remained 
non-linear,  b)  unloading  below  a  stress  level  of  -5  MPa  resulted  in 
zero  residual  strains  (i.e.  no  inelastic  response),  c)  initial  Young's 
moduli  (for  stress-strain  curve  up  to  ~2.5  MPa)  for  190  and  100  pm 
thick  specimens  are  71.1  and  69.0  GPa  and  d)  the  stress-strain 
response  of  150  pm  thick  specimen  (not  shown)  lies  in  between  190 
and  100  pm  thick  specimens.  At  the  same  time,  all  unloading 
beyond  a  stress  level  of  -5  MPa  led  to  non-zero,  positive,  residual 
inelastic  strains  and  increased  residual  strains  with  an  increase  in 
maximum  stress.  An  anelastic  response  (recoverable  strain)  is  due 
to  opening  and  closure  of  the  thermally  induced  microcracks  and 
frictional  sliding  of  irregular  crack  faces  [7  .  The  inelastic  response 
(irrecoverable  strain),  however,  is  an  indicator  of  the  permanent 
damage  in  the  material  and  could  arise  from  the  extension  of 
thermally  induced  microcracks  and/or  newly  generated  micro¬ 
cracks  during  loading.  The  microcracks  that  cause  permanent 
deformation  are  referred  to  as  mechanically  induced  microcracks. 
Based  on  the  above  observations  and  a  recent  study  on  similar 
materials  by  Pandey  et  al.,  [7],  it  can  be  concluded  that  loading 
beyond  4-5  MPa  would  introduce  irrecoverable  damage  in  the 
MMA  substrate  material  and  this  damage  is  responsible  for  a 
further  increase  in  non-linearity  and  decrease  in  the  slope  of  the 
stress-strain  curve. 

3.3.  Uniaxial  tensile  response:  multilayer  ceramic  specimens 

Table  2  presents  an  overview  of  existing  experimental  tech¬ 
niques,  which  have  been  applied  to  study  the  elastic  properties  of 
multilayer  (composites,  coatings,  thin  films  etc.)  specimens.  A  re¬ 
view  of  the  literature  suggests  that  there  are  few  established  tech¬ 
niques  available  to  measure  the  in-situ  elastic  properties  of 
multilayer  (more  than  3  layers)  ceramic  materials.  At  the  same  time, 
existing  techniques  for  testing  multilayer  specimens  are  not  readily 
applicable  for  specimens  with  porous  ceramic  layers.  For  example, 
for  porous  ceramic  specimens  the  widely  used  beam  bending 
experiment  has  limitations,  such  as  crushing  at  contact  points.  The 
corresponding  theory  for  flexural  testing  of  specimens  with  more 
than  three  layers  could  be  complex  [34].  At  the  same  time,  for  porous 
materials,  the  flexure  strength  is  found  to  be  in  between  the  tensile 
and  compressive  strength  with  tensile  strength  being  the  lowest 
[35].  A  relevant  technique  for  Young's  modulus  measurement  of 
individual  layers  in  multilayer  specimens  is  based  on  measurement 
of  curvature  change  with  deposition  of  multilayers  [48,49  .  The 


curvature  technique,  however,  would  also  require  polishing / 
grinding  of  the  thick  substrate,  since  the  sensitivity  of  displacement 
to  deposition  stress  is  directly  proportional  to  the  substrate  length 
and  inversely  proportional  to  square  of  substrate  thickness  [45].  The 
feasibility  of  this  technique  was  verified  by  using  the  Stoney  formula 
[48  ]  on  the  bilayer  beam  (Fig.  4a).  It  was  found  that  a  20  MPa  residual 
stress  (which  is  close  to  the  tensile  failure  stress)  with  a  5  mm  thick 
substrate  results  in  an  output  beam  radius  of  -680  m  and  the  cor¬ 
responding  curvature  measurement  would  be  beyond  the  mea¬ 
surement  capability  of  available  instrumentation  [46].  A  recent 
paper  by  Araki  et  al.  [47]  have  used  a  combination  of  finite  element 
(FE)  and  RUS  technique  to  determine  the  elastic  properties  of  the 
electrolyte.  However,  the  use  of  this  technique  is  limited  by  the 
number  of  layers  in  the  multilayer  specimen  and  the  corresponding 
analysis  gets  complicated  for  more  than  3  layers. 

Uniaxial  tensile  testing  tests  were  performed  on  the  multilayer 
ceramic  specimens  with  the  layers  parallel  to  the  loading  axis.  This 
loading  scheme  qualified  the  use  of  equistrain  condition  [36,50]  to 
calculate  the  elastic  modulus  of  each  individual  layer  that  forms  the 
composite.  Fig.  4a-f  shows  the  uniaxial  tensile  stress-strain  curves 
of  multiple  ceramic  layers,  starting  from  a  simpler  bilayer  specimen 
(Fig.  4a)  to  multilayer  specimens  (Fig.  4b-f).  It  is  mentioned  here 
that  for  each  multilayer  ceramic  specimen,  strains  on  all  tensile 
specimens  were  measured  in  the  gage  section  on  both  planes  (i.e. 
both  sides,  the  MMA  substrate  and  the  outermost  functional  layer) 
in  at  least  two  separate  tests.  This  is  illustrated  in  Fig.  4b-f,  where 
each  plot  shows  at  least  two  sets  of  stress-strain  curves  for  each 
type  of  multilayer  specimen.  In  Fig.  4,  the  solid  symbols  and 
markers  are  chosen  for  corresponding  stress-strain  curves,  where 
imaging  and  in-plane  strain  calculations  were  performed  on  the 
MMA  substrate  and  the  outermost  layer,  respectively.  Based  on  the 
stress-strain  response  presented  in  Figs.  3  and  4,  the  following 
observations  can  be  made: 

1.  Similar  to  the  MMA  substrate  alone  (Fig.  3),  the  tensile  stress- 
strain  response  of  bi-  and  multilayer  ceramic  specimens  are 
non-linear  (Fig.  4a-f). 

2.  The  average  failure  stress  and  strain  of  multilayer  ceramic 
specimens  (Fig.  4a-f)  is  higher  when  compared  to  corre¬ 
sponding  values  for  the  MMA  substrate  specimen  (erf  -  10  MPa 
and  ef  -  0.03%). 

3.  Failure  stresses  for  specimens  in  Fig.  4a-c  are  lower  (erf  -  20  MPa 
and  ef  -  0.06%)  as  compared  to  specimens  in  Fig.  4d-f 
(erf  -  30  MPa  and  ef  -  0.06%),  that  is,  the  specimens  that  included 
the  dense  electrolyte  layer. 

4.  For  multilayer  ceramic  specimens,  the  stress-strain  curves 
(Fig.  4a-f)  of  the  MMA  substrate  and  outermost  layer  overlap  at 
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Fig.  4.  a:  Uniaxial  monotonic  tensile  stress-strain  response  of  bilayer  specimen  i.e.  MMA  substrate  (S)  +  porous  anode  barrier  (PAB).  b:  Uniaxial  monotonic  tensile  stress-strain 
response  of  trilayer  specimen  i.e.  MMA  substrate  (S)  +  porous  anode  barrier  (PAB)  +  anode  current  collector  (ACC),  c:  Uniaxial  monotonic  tensile  stress-strain  response  of 
multilayer  specimen  i.e.  MMA  substrate  (S)  +  porous  anode  barrier  (PAB)  +  anode  current  collector  (ACC)  +  anode  (A),  d:  Uniaxial  monotonic  tensile  stress-strain  response  of 
multilayer  specimen  i.e.  MMA  substrate  (S)  +  porous  anode  barrier  (PAB)  +  anode  current  collector  (ACC)  +  anode  (A)  +  electrolyte  (E).  e:  Uniaxial  monotonic  tensile  stress-strain 
response  of  multilayer  specimen  i.e.  MMA  substrate  (S)  +  porous  anode  barrier  (PAB)  +  anode  current  collector  (ACC)  +  anode  (A)  +  electrolyte  (E)  +  cathode  (C).  f:  Uniaxial 
monotonic  tensile  stress-strain  response  of  multilayer  specimen  i.e.  MMA  substrate  (S)  +  porous  anode  barrier  (PAB)  +  anode  current  collector  (ACC)  +  anode  (A)  +  electrolyte 
(E)  +  cathode  (C)  +  cathode  current  collector  (CCC). 


lower  values  of  stresses.  The  two  stress-strain  curves,  however, 
start  to  deviate  for  stresses  beyond  3-5  MPa. 

5.  When  the  applied  stress  exceeds  ~3  MPa,  the  in-plane  strain  on 
MMA  substrate  is  lower  than  the  corresponding  strain  in  the 
outermost  layer  for  the  same  value  of  stress,  and  the  difference 
between  the  two  strain  values  increases  with  loading.  This 
result  is  consistent  for  Fig.  4a-f. 

The  above  experimental  observations  suggest  that  the  equistrain 
assumption  is  applicable  only  for  stresses  lower  than  ~3  MPa.  In 
other  words,  the  strain  in  all  layers  is  the  same  for  applied  stresses 
less  than  3  MPa.  The  invalidity  of  the  equistrain  assumption  at 
higher  stresses  could  be  explained  based  on  the  results  discussed  in 
Section  3.2,  where  initiation  of  mechanically  induced  microcracks 
(or  residual  inelastic  strains)  started  at  3  MPa  and  further  loading 
led  to  permanent  damage  in  the  MMA  substrate  and  most  likely 
propagated  into  other  layers.  The  permanent  damage  could  result 
in  the  initiation  of  strain  incompatibility  between  the  substrate  and 
additional  layers,  thereby  limiting  the  validity  of  the  equistrain 
assumption.  It  is  known  that  [50]  an  equistrain  assumption  does  not 


hold  true  if  there  is  interfacial  sliding  between  layers  and  the  latter 
can  be  enhanced  by  mechanically  induced  cracking  in  the  substrate. 
In  addition,  strength  of  multilayer  ceramic  specimens  could  be 
affected  by  residual  stresses  in  the  layers  [51,52]  and  therefore, 
addition  of  extra  layers  could  enhance  the  failure  stresses  as 
compared  to  the  strength  of  MMA  substrate  specimens. 

With  the  above  discussion,  a  new  experimental  methodology 
is  proposed  which  was  applied  to  measure  the  in-situ  Young's 
modulus  of  thin  ceramic  layers  in  the  SOFC.  This  methodology  can 
be  used  for  multilayer  systems  where  a  series  of  specimens  could 
be  fabricated,  as  mentioned  in  Section  2.  The  Young's  modulus  of 
each  ceramic  layer  is  calculated  from  tensile  stress-strain  response 
of  the  MMA  substrate  and  multilayer  ceramic  specimens  and  is 
discussed  in  the  next  section. 

3.4.  Estimate  of  in-situ  Young's  modulus 

In  the  present  investigation,  the  Young's  modulus  was  deter¬ 
mined  using  a  20  microstrain  (pe)  deviation  from  linearity  criterion 
in  the  stress-strain  curve.  The  20  ps  deviation  from  linearity 
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definition  (or  the  initial  tangent  modulus)  was  used  based  on  the 
following  considerations: 

1 )  A  20  ge  strain  deviation  from  linearity  criterion  allowed  the  use 
of  the  stress-strain  curve  below  a  stress  value  of  ~3  MPa  in  all 
cases  thereby,  ensuring  the  validity  of  equistrain  criterion.  This 
criterion  can,  therefore,  be  applied  to  determine  the  tensile 
elastic  response  of  the  MMA  substrate  (Fig.  3a)  and  multilayer 
ceramic  specimens  (Fig.  4a-f). 

2)  In  the  case  of  the  MMA  substrate,  this  definition  yields  a  Young's 
modulus  value  of  70  ±  1  GPa  (Fig.  3a),  which  was  close  to  the  value 
obtained  using  the  resonant  ultrasound  spectroscopy  (RUS) 
technique  (72  ±  2  GPa).  Values  ofYoung's  modulus  obtained  from  a 
static  technique  (uniaxial  tension)  and  dynamic  (RUS)  technique 
are  similar  for  lower  loads.  A  similar  observation  has  been  reported 
for  other  porous  and  microcracked  ceramic  materials 
[7,9,40,53,54]. 

The  average  Young's  modulus  of  composite  multilayer  specimens 
was  determined  from  all  corresponding  stress-strain  curves 
(Fig.  4a-f)  as  listed  in  Table  1  according  to  the  20  pe  deviation  from 
linearity  definition.  The  equistrain  criterion  was  later  used  on  the 
stress-strain  curve  below  ~3  MPa  to  determine  the  Young's  modulus 
of  individual  layers  according  to  the  procedure  outlined  below. 

•  For  a  190  pm  thick  bilayer  (MMA  substrate  and  PAB)  specimen, 
the  average  thicknesses  of  MMA  substrate  and  PAB  layer  were 
165  pm  (ti)  and  25  pm  (t2),  respectively.  The  20  pe  deviation 
from  linearity  definition,  in  the  stress-strain  curve,  (Fig.  4a)  led 
to  a  Young's  modulus  of  the  multilayer  composite  (with  thick¬ 
ness,  t  =  t\  +  t2  =  190  pm)  of  69.5  ±  2.0  GPa.  Now,  applying  the 
equistrain  criterion,  the  Young's  modulus  of  the  PAB  material 
(Epab)  can  be  calculated  as: 

^bilayer  =  (^1  A) ^Substrate  "P  (P^/O^PAB  =  69.5  GPa 
Epab  —  66.2  GPa 

•  This  procedure  was  repeated  for  each  of  the  other  multilayer 
specimens.  Finally,  for  a  190  pm  thick  seven  layer  specimen,  the 
average  thickness  of  MMA  Substrate,  PAB,  ACC,  Anode,  Electro¬ 
lyte,  Cathode  and  CCC  layers  were  43  pm,  25  pm,  15  pm,  7  pm, 
10  pm,  20  pm  and  70  pm,  respectively.  The  20  pe  deviation  from 
linearity  definition  to  the  stress-strain  response  (Fig.  4f)  led  to  a 
composite  modulus  of  105.2  ±  2  GPa. 

^multilayer  =  (^1  A) ^Substrate  ~P  (^2/O^PAB  T-  •  • -(^6 /O^Cathode 

+  (h/t^Eccc 

Eccc  =  127.5  GPa 


It  should  be  noted  here  that  in  each  step  of  the  calculation,  an 
average  modulus  of  the  composite  specimen  ( Table  1 )  was  used  to 
calculate  the  average  in-situ  Young's  moduli  of  functional  layers 
( "able  3).  The  use  of  outlier  values  for  composite  modulus  would 
lead  to  additional  error  propagation  in  every  step  with  cumulative 
error  37  .  An  approach  for  estimation  of  the  error  in  the  calculation 
of  the  Young's  modulus  for  each  layer  in  the  multilayer  along  with 
the  associated  assumptions  is  presented  in  Appendix  A. 

Table  3  includes  the  calculated  in-situ  Young's  modulus  (with 
the  error  estimate  based  on  analysis  in  Appendix  A)  using  micro- 
tensile  testing.  The  Young's  modulus  determined  from  microtensile 
testing  is  compared  to  the  corresponding  values  obtained  by 
application  of  the  nanoindentation  technique  on  multilayer  spec¬ 
imens  and  the  RUS  technique  on  bulk  specimens  made  with  a 
porosity  and  composition  similar  to  the  layers  themselves.  These 
values  suggest  that  the  Young's  modulus  values  obtained  using  the 
proposed  experimental  procedure  were  in  all  cases  (except  for  the 
anode)  close  to  the  corresponding  Young's  modulus  values  ob¬ 
tained  using  the  RUS  technique.  The  Young's  modulus  values  ob¬ 
tained  using  microtensile  testing  were,  however,  found  to  be 
different  from  the  corresponding  values  obtained  by  nano¬ 
indentation  on  the  various  layers.  For  example,  in  the  case  of  MMA 
substrate,  nanoindentation  yielded  a  much  higher  value  of 
250  ±  0.2  GPa  as  compared  to  70  ±  1.0  GPa  with  the  microtensile 
testing  technique.  Such  high  value  ofYoung's  modulus,  in  the  case 
of  nanoindentation,  is  obtained  due  to  indentation  on  individual 
grains  of  the  spinel  or  periclase  phase  and  was  confirmed  from 
optical  micrographs  (Fig.  5).  The  indentation  deformation  field  can 
also  interact  with  adjacent  layers  leading  to  additional  errors.  For 
example,  the  average  Young's  modulus  of  electrolyte  material 
(6ScSZ)  from  nanoindentation  is  69.8  GPa  a  much  lower  value  than 
the  225.8  GPa  obtained  by  microtensile  testing,  and  the  latter  value 
is  closer  to  the  expected  Young's  modulus  for  this  material.  It  could 
therefore,  be  concluded  that  in  the  case  of  porous  and/or  micro- 
cracked  ceramic  materials,  the  Young's  modulus  values  obtained 
using  the  nanoindentation  technique  are  affected  by  the  local 
microstructural  features  in  the  material  i.e.  grains,  grain  bound¬ 
aries,  microcracks,  porosity  etc.,  and  therefore,  should  be  inter¬ 
preted  with  caution.  The  reasons  for  discrepancy  in  the  Young's 
modulus  values  measured  by  tensile  testing  and  RUS  methodolo¬ 
gies  have  been  discussed  in  detail  in  a  recent  paper  [7  .  For  a 
multilayer  specimen,  however,  there  are  additional  variables  and 
the  reason  for  the  discrepancy  in  the  values  of  the  Young's 
modulus  measured  by  microtensile  and  RUS  testing  are  likely 
related  to  the  porosity  and  interfacial  properties  of  individual 
layers.  For  example,  there  could  be  severe  differences  in  the  defect 
distribution  (e.g.  porosity)  of  the  anode  thin  layer  and  the  corre¬ 
sponding  bulk  RUS  specimen.  It  is,  however,  noted  that  the  Young's 
modulus  value  of  the  unreduced  anode  determined  by  micro¬ 
tensile  testing  is  closer  to  the  value  reported  in  the  literature  [55]. 
The  procedure  outlined  in  this  manuscript  could  be  extended  for 
in-operando  experiments  at  elevated  temperature  and  could  be 


Table  3 

Comparison  of  in-situ  Young's  modulus  obtained  using  microtensile  testing,  nanoindentation  on  thin  layers  and  the  RUS  technique  on  the  corresponding  bulk  specimens. 


Material  (layer) 

E  (GPa)  (tensile  curves,  using  20  pe  deviation 
from  linearity)  (E,-±AEi) 

E  (GPa)  (nano-indentation) 

E  (GPa)  (bulk  value  using  RUS) 

Substrate  only  (S) 

70  ±  1.0 

250  ±  0.2 

72  ±2 

Porous  anode  barrier  (PAB) 

66.2  ±  1.5 

63  ±  11 

— 

Anode  current  collector  (ACC) 

95.3  ±  27.0 

— 

— 

Anode  (A)  NiO  YSZ 

64.5  ±  27.3 

66  ±  6 

193  ±  0.6 

Electrolyte  (E) 

225.8  ±  0.04 

69  ±  11 

208  ±  0.6 

Cathode  (C) 

112.7  ±  8.6 

43  ±  3 

143  ±  2 

Cathode  current  collector  (CCC) 

127.5  ±  0.7 

43  ±  2 

129  ±  2.0 
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Fig.  5.  Optical  micrograph  of  typical  nano  indentation  imprints  on  various  layers  in  a  substrate  supported  solid  oxide  fuel  cell  (SOFC). 


used  to  calculate  the  residual  stresses  formed  in  the  individual 
layers  during  cooling  from  the  processing  temperature  of  the  SOFC 
stack.  Finally,  the  outlined  methodology  can  also  be  applied  for 
other  functional  multilayer  systems. 

4.  Conclusions 


service  by  trade  name,  trademark,  manufacturer,  or  otherwise  does 
not  necessarily  constitute  or  imply  its  endorsement,  recommen¬ 
dation,  or  favoring  by  the  United  States  Government  or  any  agency 
thereof.  The  views  and  opinions  of  the  authors  expressed  herein  do 
not  necessarily  state  or  reflect  those  of  the  United  States  Govern¬ 
ment,  any  agency  thereof,  or  LG  Fuel  Cell  Systems  Inc. 


An  improved  experimental  methodology  was  implemented  to 
measure  the  in-situ  Young's  modulus  of  the  substrate  and  constit¬ 
uent  ceramic  layers  of  a  substrate  supported  SOFC.  The  methodol¬ 
ogy  required  the  fabrication  of  a  series  of  multilayer  ceramic 
specimens.  In  addition,  microtensile  testing  was  performed  on 
specimens  fabricated  from  the  MMA  substrate  and  additional 
screen  printed  multilayer  specimens  that  constitute  a  solid  oxide 
fuel  cell.  The  measured  in-plane  strain  on  both  sides  of  multilayer 
specimens  provided  a  means  to  ensure  the  validity  of  the  equistrain 
criterion  that  was  applied  to  calculate  the  Young's  modulus  of  in¬ 
dividual  layers.  The  values  of  Young's  modulus  obtained  from  the 
proposed  new  methodology  were  close  to  the  corresponding  values 
obtained  on  bulk  specimens  using  the  resonant  ultrasound  spec¬ 
troscopy  (RUS)  technique  but  different  from  the  Young's  modulus 
values  determined  on  similar  specimens  by  the  nanoindentation 
technique.  In  addition  to  the  Young's  modulus  value  of  the  indi¬ 
vidual  layers,  a  simplified  procedure  to  determine  the  error  in 
determined  Young's  modulus  of  the  individual  layers  with  the 
above  methodology  was  presented. 
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Appendix  1 

In  summation  notation,  the  Young's  Modulus  of  the  multilayer 
can  be  written  as: 


11  t 

^multilayer  =  ^  ^  T  •  •  • 
i=l 


(Al) 


where,  q  is  the  thickness  of  layer  i  and  other  terms  have  been 
introduced  earlier.  Using  error  propagation  rules,  the  following 
expression  can  be  derived  for  the  error  in  the  Young's  modulus  of 
the  multilayer: 


kE 


multilayer 


(A2) 


In  this  expression,  if  all  the  error  terms  are  known  or  can  be 
estimated  after  every  multilayer  test,  the  AEmuitiiayer  can  be  esti¬ 
mated  for  the  individual  tests  as  well.  Unfortunately,  the  errors  in  q, 
t  and  Ei  are  interdependent  and  estimating  them  individually  leads 
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to  an  overestimate  in  the  total  error.  In  order  to  estimate  the  total 
error  in  the  multilayer,  a  probabilistic  Monte  Carlo  simulation  [56] 
needs  to  be  performed  which  is  beyond  the  scope  of  this  manu¬ 
script.  A  simplified  procedure  is  presented  below  to  estimate  the 
error  in  the  calculation  of  the  modulus  values.  In  this  procedure,  the 
multilayer  is  imagined  as  a  bilayer  consisting  of  layer  T  and  another 
layer  comprising  layers  1  to  (i  -  1).  With  this  simplification, 
Equation  (A2)  can  be  rewritten  as: 


multilayer 


(A3) 


where,  the  terms  with  subscript  (1  ,i  -  1)  corresponds  to  properties 
of  all  but  the  ith  layer  and  the  latter  is  denoted  by  subscript  (i).  For 
calculations  associated  with  the  layer  consisting  of  (1 ,  i  -  1)  layers, 
it  was  assumed  that  error  in  modulus  measurement  was  2.5%  and 
for  all  conditions  the  total  error  in  thickness  measurement  was 
assumed  to  be  0.4%,  i.e.  +  A t/t  =  +  A t/t  =  0.004 

and  A E^i-i /£!,!_!  =  0.025.  Assuming  higher  errors  would  lead  to  a 
negative  AEmuitiiayer  which  will  be  physically  unrealistic.  Since  all 
the  terms  are  known  in  Equation  (A3),  inputting  the  experimentally 
determined  error  in  the  modulus  of  the  multilayer,  which  is  the 
error  in  the  modulus  value  of  the  composite  in  Table  1,  will  give  the 
error  in  the  calculated  modulus  of  the  individual  layer,  A £/.  This 
procedure  was  applied  for  all  the  layers  and  in  this  manner,  the 
Ej±AEj  values  are  reported  in  Table  3. 
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